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Abstract 

In connection with the contribution "Quantum Condensates in Nuclear Mat- 
ter" some problems are given to become more familiar with the techniques of 
many-particle physics. 

o 

; 1 Bogoliubov transformation for superfluid state 

S3 . 

The entropy operator S, which is related to the (grand canonical) statistical 
operator via p = cxp[— S/fcs], has in second quantization the form S = So + 
S1+S2 • • • where Sq is a normalization constant, the single-particle contribution 
has the general form 

Si = 22s(i,j)a}aj +2_ d d{i,j)a i a j + d*(i, j)a]a\i 

T I ij ij ij 

where the Lagrange multipliers s(i,j),d(i,j),d*(i,j) are determined by the 
given averages when maximizing the entropy, the two-particle contribution and 
higher terms are neglected in the mean-held approximation considered here. 
Assuming that Si is hermitean and that the normal term s(i,j) is diagonalized, 
we have 

CO ■ p p p 

if we furthermore assume that the condensate is formed with a given spin and 
momentum state for the two-particle system (p denotes momentum and spin 
. quantum number). In particular we can take s(p) = p 2 /2m a — fx a , and the pair 

I amplitude couples p = p,a with the state p = p, —a 

Problem 1: 

•rH , 

' Find the transformation which diagonalizes S\ . 

Solution: 

The Bogoliubov transformation for the creation and annihilation operators reads 

a p = u p b p - v p b' p a\ = u* p b\ - v*b p , (1) 
a p = Upbp + v p b p a P = u* p bp + v*b p , (2) 

where the new operators obey the anticommutator relations 

{b p ,bl,}=S p y , {b p ,b p ,} = {bl,bl,} = (3) 
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In order to get a canonical transformation, i.e. the anticommutator remains 
unchanged, we have to claim {a p , a p ,} — 6 PjP /. 



{a p ,a f p ,} 



m 



Upb p u*,b p , - Upb p v*,bp, - VpbfeUpty + v p b' p v*,bp' 
+u*,b p ,u p b p - v*,bp,u p b p - u*,b p ,v p bp + v p ,b p >v p bp 

Op| 2 + \v p \ 2 )5 p , p/ 



So it follows 



\ Up \ 2 + \vp\ 2 = l . (4) 

We now list the prefactors in Si where we make use of the anticommutator 
relations ([3]). In order to achieve a diagonal transformation the inner columns 





bpbp 


bpbp 


b f p bl 


bpbl 


s(p) 

dip) 


CLpClp 


~\v P \ 2 

UpVp 


-Kl 2 

UpVp 


-u p v* 
up 


-VpU* 

-4 


d*(p) 

sip) 


CLpClp 
dp dp 


u*v* 

Kl 2 


-u*v* 

K\ 2 


~v* 2 
p 

-u p v* 


u*p 2 

-UpVp 



Table 1: Coefficient matrix of the canonical (Bogoliubov) transformation 

of Table Q] have to vanish, i.e. after substituting s(i) = e p and dip) — A p we 
have 

- 2e P iupV* p + u* p v p ) + A P iu 2 p - v 2 ) + A* p iu* p 2 - v* 2 ) = (5) 
One can easily check that 



1 


1 + 


e P 




72% 


y/4 + 


A 2 

p 


1 


1 - 


e P 




V2% 


^1 + 


A 2 

p 



(6) 
(7) 



solve the equations Q and ([S]) with e p , A p being real functions. Furthermore 
we obtain 

Slip) = spl + ^l (blbp + blbp- 
which is now diagonalized. 



2 Transition from BCS to BEC 
Problem 2: 

Give the form of the wave function for the two-particle wave function including 
Pauli blocking for the Yamaguchi interaction (isospin singlet) at E = 2/z. How 
changes the wave function in coordinate-space representation if we cross over 
from low densities (deuteron) to high densities (Cooper pairs)? 
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Solution: 



The Schrodinger equation for the two-particle problem including Pauli blocking 
reads 

[ j E(1)+^(2)-£„p]^p(12)+^[1-/(1 , )-/(2')]V(12, l'2')Vnp(l'2') = 0, (8) 

1'2' 

where E(l) = E(pi) — p\/2m\ denotes the kinetic energy of the single-particle 
state {1} = {pi, 01, Ti} abbreviating linear momentum, spin and isospin, respec- 
tively. The Fermi distribution function /(l) = {exp{[E(pi) — P i]/T} + is 
characterized by the temperature T and the chemical potential 
The Yamaguchi potential is 

y(12,l'2') = -jjW(( Pl - P2 )/2)W(( Pl ,-p 2 ,)/2)5 Pl+P2 , Pl , +P2 , 

with the formfactor W(p) = (p 2 + 7 2 ) -1 , the effective range 7 = 285.8484 MeV 
and the effective coupling constant Xs=i,t=o = 0.4144 fm -3 for the spin triplet, 
isospin singlet channel and M is taken as the average nucleoli mass. 
After restriction to a two-nucleon system at rest in an isospin symmetric medium 
(/ii = fi2 = H\ ™>i = m2 = m )i the wave function satisfying the Schrodinger 
equation (jHJ) reads 

^(p) = c(r >/ i) p2/m ^ (I>) 0) 

with the normalization factor c„ (T, /i) . Inserting © into © one obtains the 
implicit equation for the binding energy E a (T, //) 

\ aSr W\p')[l-2f{p')] 

m Z-f p iy m - E a (T,fi) ' 1 ' 

It is evident that the Pauli blocking factor in the kernel of (fTU)) generates the 
T- and fi- dependence of the binding energy. 

The wave function in coordinate space is obtained by Fourier transformation, 

1 f 00 

= — — / dpp Vv(p)sin(p r), (11) 
^ r Jo 

with the final result 



C m e -V-mE a r _ e ~ir 

Mr) = 4xtf + mE 9 ) r ' (12) 

A zero binding energy defines the transition from a bound state to a scattering 
state (Cooper pair). The corresponding transition from negative to positive 
energy eigenvalues entails a character change in the wave function to oscillatory 
behaviour. 
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The Schrodinger equation for the scattering problem (with the same restrictions 
as above) reads pQ 

(k 2 - P 2 + ie)i>(p) = -KW( P ) ]T[1 - 2f{p'))W{p')^{p') , (13) 

p' 



with the solution 



fc 2 + 7 2 

2ir 2 " v ' 7 2 + p 2 k 2 — p 2 



= S( P k) - -£J- ■ F a (k) ■ —n——n ■ TT—T— (14) 



where !Fa{k) is the scattering amplitude and can be obtained by inserting (|14[) 
into ([lU)) . When we neglect Pauli-Blocking, T a (k) takes the form 

Including Pauli-Blocking, the numerical solution is given in [2J. 

Equivalently, one can reprsent the solution in terms of the phase shift <5 CT (/i,T), 

given by 

jr (fc) = Sm ^ = 1 Qg-j 

A: — zfc + A; cot ho- 

using (IT5t and _E = k 2 /m, one obtains the explicit form for the scattering phase 
shift 8 (T {E) from 

cot5 a {E) = - 7 ± r {l-x-y{l + xf) (17) 

where we have introduced dimensionless variables y = 7 3 /(7r 2 A CT ) and x — 
E a /E with E — 7 2 /m. 

In order to discuss qualitatively the Mott-Effect of vanishing bound states we 
consider (ITU)) without Pauli-Blocking but with variable potential strength. Ac- 
cording to the Levinson theorem, at the critical coupling strength for the dis- 
solution of the bound state, the scattering phase shift has to jump by it at the 
threshold. This critical coupling strength can be calculated by solving ([TU]) or 
(JTTJ) respectively. 

Let us consider the case when E tends to zero (x —> 0). (IT71) gives 

lim cot^ = 77^(1 - y) U-+o • (18) 

Using (I10[) and setting the binding energy to zero, we obtain the critical coupling 
A c = 7 3 /7r 2 which corresponds to y c = 1. When approaching this critical 
coupling from both sides, we get 

lim cot SJE) =±oo (19) 

A->A c ±0 

which corresponds to a jump of the scattering phase shift from it to 0, in accor- 
dance with the Levinson theorem (illustrated in Figure [1]). 

We now want to consider the Pauli-Blocking energy shift by the example of 
the Deuteron. 
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Figure 1: The scattering phase shift as a function of energy for several values 
of the coupling. Ac is defined in the text. 

3 Pauli blocking for Gaussian bound states 
Problem 3: 

Construct a wave function for the bound states of 2, 3, 4 nucleons by minimizing 
the energy (Yamaguchi interaction) with respect to Gaussian wave functions! 
Calculate the Pauli blocking energy shift within first order perturbation theory! 
What is its dependence on temperature? 

Solution: 

The four-particle Schrodingcr equation reads 



£ nP V„p(1234) = [E{1) + E{2) + E{2>) + E{A)]^ nP {l2M) 
+ E {[l-/(l)-/(2)]ni2l'2%3'& 



'44' 



1'2'3'4' 




3.1 a-like clusters 

We make a Gaussian ansatz for the wave function, i.e. 

i>(n,r 2 ,r 3 ,u) = N^e-^-^Mr.-KfHr.-nfHr.-nf)^ 
where R = (ri + r 2 + r 3 + Ta)/A, or in Jacobi-coordinates 

where we have used the transformation rules 

1 k n n—1 

i—i i—i i—i 



(20) 



(21) 



(22) 



The Fourier transformation of the wave function reads 

iP(h,k 2 ,h,k 4 ) = J V(ei,6,6,i?)e- t( « l9l+ « 292+?393+Kii) d 3 arf 3 6^3d 3 i? 

= 7V X (X) e -^(T^+| 92 2 +l^)^ fci+fe2+fc3+fe4 (23) 
where K = k\ + k 2 + k 3 + fc 4 and qi are the conjugate momenta to 

3.2 Normalization 

The normalization factor can be determined as usual by calculating the expec- 
tation value of the square of the wave function: 
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7V2 



fin 



(2tt) 3 

a 



e~^i'ii+iii+iii)d 3 q 1 d 3 q 2 d 3 q 3 

3 



(27T) 

fi 



2_ J e-^^+h 2 2 + U) dqidq2dq3 

/oo 
-oo 



(2tt) 3 V a 2 V 2 3 4 

fip 

(2^) E 



3 v^V 



so it results 



(2tt) 3 



fin 



(24) 



3.3 Kinetic energy term 

(E kln ) = N 2 ^{ k l + k l + k l + k l) \^{ki,k 2 ,k 3 ,k 4 ) 



k 1 ,k 2 ,k 3 ,k4 



6 



/Ti C/n Ho V ' 



d 



e 2 



■(!9?+|92 2 +l« 3 2 ) 



2™<9(-±a 2 ) 
km 2m a da N 2 

pK , J_A 

where we explicitly separated the motion of the total momentum K. 
3.4 Potential energy term 

For a separable potential we can calculate the potential part of the energy as 



(Epot) — -3- 



n 

, A S + A t 



w(q 1 )w(q[)ip(q 1 ,q 2 , q 3 , K)ip(q[, q 2 , q 3 , K) 



91,92 , 93, K,q[ 



]T w(qiMq[)N 2 | x(K) | : 



xe 



= -3 



91,92, 93, K,q[ 



A, + A t 1 



1 J2 w (lMl')e-^ 2+ "' 2 



9,9' 



L 9 



As an example we assume the potential to be of Gaussian type w(q) = exp{— q 2 /j 2 }, 
this leads to 



(E po t) — -3 



A s + At 



(2tt)3 e -#V 



^ 9 d 3 



-24 



Ag + A t a 3 y/i 
(2»)« "'" ~" 



3 ' 



So finally the total energy is given by 

9 1 

2 ma 2 



(■Hint) - (H) - £^„ 



24 



As + A t a 3 ^7r 



(25) 
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3.5 Model calculation 



We can now minimize this expression with respect to the variational 
parameter a and fit the binding energy for the a particle by adjusting our 
model parameters. We obtain 

A t = 1317.8 MeV fm~ 3 , A s = 667 MeV im~ 3 

and therefore 

a = 8.288 ■ 1CT 3 MeV" 2 , 

which gives us a binding energy of the a particle of —28.2 MeV. 

For the other bound states with less nucleons, this works in a similar way. 

3.6 Pauli blocking shift 

Now we want to calculate the Pauli blocking shift due to finite temperature. 
We will consider a bound state of 2 nucleons. For 3 and 4 nucleons this works 
similar but is more cumbersome. From (|10[) we know that the exact wave 
function in the vacuum satisfies 

\ a v w 2 ( P ') 

m^p' 2 /m-E% [ ' 

pi 

In perturbation theory we expand the binding energy around its vacuum 
value, i.e. \E a \ = E® — AE Pauh (T, //) where the minus sign is because the 
binding energy is negative (E a = — \E a \). 

1 _ ^ W 2 (p')[l-2f{p')] 



p 



f p' 2 /m + E°- A£; Pauli (T, n) 



Ky W 2 (p>) 1-2/0/) 

m ' J 



p' y ' p>y m +E« 



. ^ E J^( 1 _ 2/M+(1 _ yW) ^!e>) +0((aE 

m p' z /m + E% \ p' z /m + E% v 

Dropping higher order terms we arrive at 

Ap Paul irT lA f\m\ 2 (P 2 /rn + E°)2f(p)p 2 dp 

l\E = , (27) 

J \ip{P)\ P dP 

we have assumed 2f(p') to be small against 1 and used the expression for the 
wave function ip(p) given in Eq. ([9]). 
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